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for comparing the usual abelian Higgs approximation used in cosmology to realistic situations. 
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1 Introduction 

Decade long improvements in experimental data led to the conclusion that only those Grand 
Unified Theories (GUT) involving some amount of Supersymmetry (SUSY) were acceptable [1-9]. 
The vacuum structure of these theories implies that they should have produced topological defects 
during their successive steps of Spontaneous Symmetry Breaking (SSB), including monopoles and 
cosmic strings [10, 11], a phase of cosmic inflation being then necessary to dilute the former. If 
we consider furthermore a E—term hybrid inflation scenario [12-18], most of the SSB schemes 
lead to the formation of cosmic strings at the end of inflation [19]. So, constraining the string 
energy per unit length, e.g. through Cosmic Microwave Background (CMB) observations [20-24] 
provides a general way to constrain GUT themselves. 

The structure of cosmic strings forming at the end of inflation has already been studied in 
details, see Refs. [25-32]. These works considered models where only the minimal field content 
necessary to form a string was introduced, and where the scale of formation of strings is the 
only dimensionful parameter. The aim of the present paper is to investigate the structure of 
cosmic strings in a realistic GUT context, i.e. starting from the complete GUT field content, and 
analyzing how the Higgs fields implementing the SSB scheme form the strings. This is done in 
F-terms models, focusing only on the bosonic part of the supermultiplets. 

Considering a general SUSY GUT, we identify the minimal field content sufficient to de¬ 
scribe the realistic string structure. This minimal structure involves all the fields which take non 
vanishing Vacuum Expectations Values (VEVs) at the end of the SSB scheme, and thus are singlet 
of the Standard Model (SM). The energy par unit length will be modified by this condensation of 
several Higgs fields in the core of the string. The additional fields also give natural candidates to 
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carry bosonic currents [33-37], and we expect their superpartner to carry zero-modes fermionic 
currents [33, 38-42]. Furthermore, the intercommutation process [43-48] and the cusps evapora¬ 
tions [49-53] can be qualitatively modified due to this extra structure. The modification of any of 
these properties can have a major impact on the cosmological consequences of cosmic strings. In 
this paper, we give a complete description of the realistic microscopic structure of cosmic strings, 
and we perform a perturbative study of the modification of their energy per unit length from 
standard toy models. This gives a first step for the study of the other phenomena mentioned 
above. 

For this purpose, we give an ansatz and boundary conditions for such a minimal structure, in 
the case of Nambu-Goto abelian strings. The conventions and normalizations chosen are such that 
if all the fields were to decouple from the string-forming Higgs, one would recover the standard 
abelian Higgs model. Two different classes of strings are discussed, only depending on how the 
Higgs helds of the GUT implement the SSB scheme, referred to as single and many-field strings. 
We perturbatively evaluate the modification of the energy per unit length from standard toy 
models due to this complex structure, taking into account the numerical factors appearing when 
one writes the complete formulation of the GUT. 

In Sec. 2, we specify the SUSY GUT studied while setting the notation, and briefly describe 
its SSB scheme in parallel with the inflationary process. In Sec. 3, using properties of the GUT 
reviewed in the first section, we present the abelian cosmic strings considered, discuss their 
minimal structure, and distinguish the two categories of strings. Finally, in Secs. 4 and 5, we 
propose ansatze and perform a perturbative study of their properties. 

2 Theoretical framework 

2.1 Field content and snperpotential 

We consider a general SUSY GUT associated with a gauge group G. Usual examples are based on 
either SO(10) or SU(6). The spontaneous symmetry breaking down to the SM will be implemented 
in the context of a F-term hybrid inflation [12-18]. The field content of the theory includes a 
set of chiral supermultiplets and a gauge supermultiplet associated with the generators of G. We 
restrict ourselves to the bosonic sector of the model, and so only write down the scalar part of 
the chiral supermultiplets. As we study a F-term theory, we assume the F-terms are identically 
zero, with no Fayet-Iliopoulos terms [54]. It will add some constraints to the fields. 

In order to implement hybrid inflation, it is a necessity to have at least two helds 'S and 
5] in complex conjugate representations, which have a coupling term with the inhaton in the 
superpotential [12-14]. The inhaton is assumed to be a chiral supermultiplet of scalar component 
S, singlet under G. In order to reach the SM symmetry and have a phase of inhation washing 
away monopoles, one needs to have more than one SSB step; this means other helds should be 
present. The other chiral supermultiplets are denoted for the helds in real representations, 
and and for the helds in complex representations. 

We assume the most general superpotential taking into account all the above chiral super¬ 
multiplets, supplemented by a specihc term for the inhaton S to implement hybrid inhation. It 
is, assuming explicit summation on all repeated indices, [54] 

W = j -h 

+ Pxy'^^x^y + Pxy'^^x^y + -^^xyz^x^y^z T kS{'S'S — M^), (2-1) 

the last term actually implementing hybrid inhation. In this equation, the label x, y and z can be 
either I, i or 7. All the coefficients which appear in addition to the helds are complex constants, 
and /3xy, (5xy and Xxyz are totally symmetric in their indices. The constants k, M, ms and 
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the diagonal elements of the mass matrices can be set real and positive by redefinition of the 
phases of the fields. Depending on the explicit choices of representations, the coefficients m/j, 

Vx, Pxy, Pxy and Xxyz will be non-zero only when they allow to build gauge singlets. We 
did not include terms like and for the sake of simplicity; their contribution to the 

macroscopic structure of the string is shortly discussed in Sec. 5.2, and found similar to that of 
terms like 

In this work, we chose to consider separately the GUT and inflation parts, being as realistic 
as possible for the former and setting up the simplest possible model for the latter. This is the 
approach found also e.g. in [55]. This simple inflation term is sufficient to reproduce the standard 
inflation phenomenology. Besides, different kinds of terms such as or 5'$^ generate mass 

or quartic terms for the inflaton in the scalar potential and can thus spoil the inflation. They 
must therefore be considered with care [55]. 

In the example of a SO(IO) GUT, the fields X) and S are often taken to transform as the 
126 and 126 representations, which are the lowest dimensional complex conjugate representations 
which are safe, i.e. permitting i?-parity conservation at low energy to ensure proton stability [56]. 
The other fields used to implement the SSB can be for example a 210 and a 10 [8], or two 54, 
two 45, and two 10 representations [57]. The whole expression of a SO(IO) model, with the 
normalizations and conventions of the present paper, is given in Ref. [58]. 

2.2 Lagrangian of the bosonic sector 

Let us write down the general form for the Lagrangian of the bosonic sector, while setting the 
notation. In the following, we take the signature of the metric to be -|-2, and label with latin 
indices a, 6, ... the generators of G, with gauge coupling constant g. The kinetic part of the 
Lagrangian is thus (with implicit summations on i, i and I) 


- - (V^5)*(V^5) - (2.2) 

with 

D^X = (V^ - igA^^T^)X, (2.3) 

being the relevant operators encoding the action of the generator labeled by a on the field 
X. From now on, we will denote by X a generic scalar held when unspecihed, i.e. X £ 
{X), S, $/, $ 4 , S}. The strength tensor is dehned in the usual way, 

+ gf\,A^^Al, (2.4) 


with /“he the structure constants of G. 

The potential term is constructed from the F-term, obtained by taking the derivative of the 
superpotential with respect to the chiral supermultiplets 


This yields the terms 


Fx = 


dW 


(2.5) 


Fx; — “1“ Pxy^x^y'i 

“ 1 “ Vx^x “ 1 “ -|- pxy^x^y} 

= mjj^j + ?7/SS + 2pjy^^y + 2Pjy'S^y + Xjyz^y^z^ 
F$^ — ‘^Piy^^y ‘^l3iy^^y Xiyz^y^z-) 

F$- TTlji^ j “h “h 2Piy^^y H" X^yz^y^ Z"! 

F 5 = ac(SI]-M 2 ), 
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Fx being in the conjugate representation of X. These T-terms are discussed in more details in 
Sec. 4.2. The scalar potential is finally obtained from these terms through 

V = (2.7) 

X X 

where Vx > 0 and V > 0, and where we use additional bold symbols for the T-terms to remind 
that they are not singlet of the gauge group in general. Note that we did not include in V the 
D-term contribution, since these terms identically vanish in a T-term scenario, and thus play 
no role in the dynamical study of the fields. However, they are indeed taken into account by 
imposing some constraints which are discussed in Sec. 4.1 and 5.1. 

Finally, the full Lagrangian density is derived from these two terms 

C = K-V, (2.8) 


i.e. by adding Eqs. (2.2) and (2.7). 


2.3 Hybrid inflation and SSB scheme 


It has been shown in [19] that the formation of cosmic strings at the end of an hybrid inflation 
phase is essentially unavoidable during the final stage of a SUSY GUT symmetry breaking. It is 
this category of strings that we discuss in the present paper. In such models, all the monopoles 
which may have been produced in a previous phase are washed out during inflation. 

The matter content we introduced previously implements the SSB scheme down to the SM 
in at least two steps : 


G 





> Gsm X Z2 


(2.9) 


In many cases, the symmetry that is broken at the last step contains but since there 

are no general constraints about it, it is better to leave it arbitrary. 

All the non vanishing VEVs after the end of inflation have to be singlet under the SM gauge 
group, otherwise the vacuum would have non vanishing quantum numbers under this symmetry. 
We also assume there is no symmetry restoration, i.e. all fields acquiring a non-zero VEV at a 
given stage keep it non vanishing at later stages. So, we can restrict the study of the SSB scheme 
to SM singlets only. 

At the onset of inflation, we can assume an initially very large value for the inflaton S in 
comparison with all the other fields, as is expected for chaotic inflation. To minimize Vs ~ 
and ~ |«:S5p, the fields S and H must take a vanishing VEV. The terms f^xy^x^y niust 
be in the same representation as H for the terms Pxy^x^yF^ to be scalars. So, if they take a 
VEV before the end of inflation, all the symmetries broken by S after inflation would already 
be broken at this current step, which is contrary to our assumptions. Thus, they cannot take a 
non-zero VEV before the end of inflation. This property and the same reasoning for Pxy^x^y 
show that Vb = 0 and Ug = 0 before the end of inflation. This also ensures that all the terms 
implying the inflaton S identically vanish in the potential at tree level. 

Einally, we can assume that at the onset of inflation, the VEVs verify that all the potential 
terms except Vs are zero; it is indeed the global minimum for the potential taking into account 
the constraint on Vs- Different field configurations give this minimal value, including that with 
all the fields having a vanishing value^. So, it gives several sets of solutions 


{(^(-)) — Oj (^{-)) — (^a;,(-))}5, (2-10) 

the index (—) meaning that we consider the set of VEVs before the end of inflation, while 5 labels 
the different sets themselves. 

^We assume that other configurations can exist, since the potential conditions give fourth order polynomial 
equations in the fields. 
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At the end of inflation, the fields reach a global minimum for the potential, meaning all 
F-terms contributions must independently vanish. Setting 1^ = 0 then implies 

(S)(S) = Af2, (2.11) 

and the fields X) and X) take non zero VEVs as expected. The vanishing conditions of all the 
other potential terms finally give several sets of solutions 

{(^(+))5 (^(+))) ‘ S '(+)) (^ a ;,(+))} 5 ', ( 2 - 12 ) 

the index (+) denoting that we consider the set of VEVs after the end of inflation, while 6' labels 
the different sets of solution. 

During inflation, the potential for the inflaton yields V = Vq + quant.corr., where Vq = 
is the value of the potential at tree-level. Due to the quantum corrections, the value of the inflaton 
field will then slowly roll, until it meets its critical value, thus ending inflation (see Ref. [55] for 
an explicit example). 

As we study strings in models which are relevant in a particle physics and cosmological point 
of view, we assume that the superpotential and field content used imply at least one appropriate 
SSB scheme. For these schemes, the set of VEVs at the end of inflation defines the SM gauge 
group, no harmful topological defects are produced, and the stability of the inflationary valley is 
ensured (see the associated discussion^ in Ref. [55]). We consider from now on such a SSB scheme 
and the associated non vanishing VEVs, those being described as in Eqs. (2.10) and (2.12) without 
considering the labels 5 and 8' anymore, i.e. explicitly assuming a specific set in each ensemble. 

Such a SSB scheme in parallel with the inflationary process in a given SO(IO) GUT is 
described in Ref. [58] . 

2.4 Description with the restricted representations 

As we work with large dimensional representations, it is useful to consider their branching rules 
to simplify the description of these fields. Indeed, we only need to work with fields behav¬ 
ing as singlets under the SM gauge group to describe the whole SSB scheme. For instance, 
the 210 representation of SO(IO) contains three such restricted representations, included in 
its (1,1,1), (1,1,15) and (1,3,15) representations of its decomposition under the Pati-Salam 
group [SU(2)ixSU(2)ijxSU(4)] [59]. The 126 and 126 representations only have one such sub¬ 
representation, contained respectively in their (1,3,10) and (1,3,10) restricted representations 
under Pati-Salam. 

Let us consider such a VEV singlet. For a given field $ 3,5 we ascribe an index a to describe 
the different sub-representations transforming trivially under the SM, and write the associated 
VEVs as 

(*x,a) = (x^) (^a:,a)o> (2-13) 

without implicit summation, and where (l>x,a is a complex function of space-time and (^a;,Q;)Q a 
constant normalized vector in representation space ((^a;,o)o(^ 3 :,o)o = !)• Using these notations, 
the complete part of the field which is singlet under the SM can be written as the combination 

= (2-14) 

a 

This procedure reduces the fields description to only a few complex functions. 

^It was also argued in Ref. [55] that one should verify that the set of VEVs before the end of inflation is close 
in field space to only one set of VEVs after the end of inflation, since otherwise the possibility that the fields take 
different sets of VEVs at the end of inflation could create domain walls. 
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3 Abelian cosmic strings 

3.1 Strings studied 

Let us focus on the strings created at the last step of SSB of the GUT, lowering the rank of the 
gauge group by one unit and ending the inflation phase. Several kinds of strings can appear, 
depending on the quotient group H ~ G'/Gsm- As the SSB lowers the rank of the group, H 
must contain at least one U(l) subgronp, and we will focus on the Nambu-Goto abelian strings 
which form at this step, associated with this abelian generator. When H is larger than U(l), non 
abelian strings conld also form [26-28]. However, and since we want to constraint all GUTs, we 
restrict attention to the minimal U(l) case, any other kind of strings tightening the constraints. 
We will denote U(l)str this particnlar snbgroup and r®*’’ the associated generator. 

These strings cannot be connected to monopoles. As shown in Ref. [60], the strings are stable 
with respect to breaking into monopoles when the scale of formation of these monopoles is higher 
than the scale of formation of the strings. It is the case here, as we assumed that only strings 
form at the SSB considered. So, these strings could connect only to pre-existing monopoles, and 
those have, by construction, already been washed away during the inflation phase. 

In what follows, we nse a set of cylindrical coordinates (r, 6, z, t) based on the location of 
the string, and taken to be locally aligned along the z-axis at r = 0. We also focus on strings 
with fields functions of r and 6 only, and consider uniquely the bosonic structure of the string. 
It means that we will not consider in this paper any currents in the core of the string. This 
possibility will however be discussed in the following section. 

3.2 Minimal structure 

We have to determine which sub-representations are sufficient to describe the structure of the 
string. On the one hand, the fields take at infinity the non vanishing VEVs defining the SM 
symmetry, so all the restricted representations non singlet under this symmetry take vanishing 
values far from the string. On the other hand, the sub representations non singlet nnder the SM 
appear at least in a quadratic form in the potential. Indeed, a potential term containing only one 
such field would be charged under the SM. Both these resnlts imply that an ansatz where all the 
sub-representations charged under the SM take an identically vanishing value is solution of the 
equations of motion with the boundary conditions at infinity. 

To understand in another way this ansatz, one can consider the static configuration which 
minimizes the potential at the center of the string. Such a configuration is given by the set of 
non vanishing VEVs singlet under the SM defining G' in the SSB scheme, see Sec. 2.3. These 
VEVs are not charged under U(l)str) which wonld otherwise already be broken at this step. 
Using arguments discussed e.g. in Ref. [33], we expect fields in the string to take intermediate 
values between the configurations which minimize the potential at the center of the string and at 
infinity, depending on the competition between kinetic and potential terms. Then, as both these 
confignrations only contain singlets of the SM, one expects an ansatz where all the other fields 
take an identically vanishing value. 

This particular ansatz, which we consider in the following, is what we define as the minimal 
structure. As the configurations which minimize the potential at the center of the string and at 
infinity are a priori different since they are solutions of two different quartic polynomial equations, 
all the fields taking non vanishing values at the last step of SSB condense in the string, and thns 
have to be taken into acconnt. Note that a complete study of the stability of such an ansatz 
shonld be considered in each given model. 

The minimal structnre for cosmic strings, where several bosonic fields condense in the string, 
should modify most of the properties of these objects. Some of these properties are described 
below. Eor instance, it becomes mandatory to describe the microscopic structure of the strings 
by taking into account the condensation of these additional Higgs fields. This then permits to 
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evaluate the modification of the energy per unit length due to this complex structure. As most 
of the observational constraints on this macroscopic parameter are rather stringent, every sizable 
modification of it in a given model could rule out this model. 

Several other properties of cosmic strings will be modified by the fact that the actual struc¬ 
ture is more involved than that of the toy models usually considered. First, the condensation of 
several Higgs fields in the core of the string gives natural candidates to carry bosonic currents [33- 
37]. Also, and as we work in a SUSY framework, we can expect the superpartners of these Higgs 
fields to build fermionic currents via their zero modes [33, 38-42]. This more complex structure 
can also qualitatively modify the intercommutation process [43-48] , which has a tremendous im¬ 
pact on the temporal evolution of the cosmological string network, and thus on the consequences 
on the CMB [10, 61-65]. Furthermore, the link made between the fields forming the string and 
the particle physics model used allows a more detailed examination of the cusps evaporation 
phenomenon [49-53]. A modification of each of these properties can have major consequences on 
the cosmological implications of cosmic strings. 


3.3 Equation of state, toy model limit 


Since we explicitly assume currentless strings, nothing in the configuration we are interested in 
can depend on the internal string worldsheet coordinates, here locally z and t. We have 

= + (3.1) 

yielding T** = = C. Then 

U = 2tt [ rdr r“ = -27r f rdr = T, (3.2) 


i.e. the Nambu-Goto equation of state, Lorentz-invariant along the worldsheet. Thus, the only 
parameter of interest is the energy per unit length U defined in Eq. (3.2). This parameter, 
frequently denoted by // in the cosmology literature, is directly constrained from, e.g., CMB 
observations: from Planck and WMAP data, there is a constraint of G^/(? < 3.2 x 10“'^ at 
95% confidence level for the abelian Higgs model [23]. Thus, every modification of U due to the 
realistic structure of the strings will have to be compared to the already stringent observational 
constraints, and could rule out the associated model. 

To translate the standard abelian Higgs model to a F-term SUSY formalism, three fields 
are necessary, S, S and S, the first two fields having opposite U(l) charges, the last one being 
uncharged. The superpotential is 

W = kS (sS - . (3.3) 


Assuming that S identically vanishes, it yields as expected the standard U(l) string model for E 
and S (see for example [26, 29, 34]), with^. 


£ = -(F^S)t(F'^S) - (F^S)t(F^S) - - 


SS - M" 


(3.4) 


We recover explicitly this limit from the realistic model when all the parameters but k and 
M go to zero, and considering only U(l)str defined in the previous part. Note that it is the case 
due to the previous normalizations and conventions, see Sec. 4 for the details. This well defined 
limit must be verified when using results from the abelian Higgs model, since important GUT 
numerical factors can appear. For this well known problem, we have S ~ S ~ M, a characteristic 
radius for S of and a characteristic energy per unit length Uq ~ [11, 26, 34, 58]. 


®Note that if one is more familiar with the model containing a single kinetic term of the form (_D^E)(Z)''E), usual 
for non supersymmetric theories, a link between both models can be easily performed. For this, it is sufficient to 
introduce new variables of the form E = %/2E, M = \/2M, and k = k/ 2. Indeed, one sees that the superpotential 
in term of the new variables remains the same, while the factor 2 of the kinetic term will disappear. This is only 
valid in a P-term scenario, where E and E are complex conjugate due to the _D-term condition (see 4.1). 
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3.4 Two classes of strings 

We assume for the entire paper that 'S and S only have one sub-representation singlet of the 
SM, charged under U(l)str- If h was not the case, we could treat the additional representations, 
charged or not under U(l)str; in the same manner than the additional fields 

We distinguish two classes of strings. For the first kind, no other fields are charged under 
U(l)gtr. We call these strings single-field strings, since only one field (and its conjugate) is directly 
coupled to the gauge field. They are discussed in Sec. 4. For the second case, other fields can be 
charged under U(l)stri and we naturally call them many-held strings, see Sec. 5. This property 
only depends on the held content of the GUT and how it permits to implement the SSB scheme. 

In the case of a single-held string, no /S-couplings, e.g. in can appear in the super¬ 

potential between the restricted representations singlet under the SM, since such terms would be 
charged under U(l)str- It is not anymore the case with a many-held string. 


4 Single-field strings 


4.1 Ansatz and boundary conditions 

According to our dehnition of single-held strings, we consider the case where only the restricted 
representations of S and X) are charged under U(l)str- We normalize their charges to gs = 1 
and qf^ = —1, the charges being dehned by identifying tx = qx^^x for an abelian generator 
in Eq. (2.3). To describe the different sub-representations singlet under the SM, we use the 
decomposition of Eq. (2.14), which gives 


(5](r,0)) = a(r,0)(S)o, (4.1) 

and 

(E(r,0)) = d(r,0)(S)o, (4.2) 

with (e:)J(e:)o = 1 and (s:)^ = (S)J. 

The H-term condition associated with i.e. [54] 

= -gY,{X^TfX) = Q, (4.3) 

ensures that a and a have the same norm. In addition, as the phase of the inhaton S have been 
rephased in order to make M real, it ensures that the global minimum of the potential is reached 
when SX) = G M. Both these results impose that S = which finally gives a = a*. 

The string itself is defined through [11, 61] 


(^)(r=0) — (^)(r=0) — 0- 


(4.4) 


Thus, we can introduce an ansatz similar to the Nielsen-Olesen solution, with integer winding 
number n : [11, 26, 27, 34] 


a = /(r)e*"®, 

4’x,a — 4^x,ai'l ")) 

S = S{r), 

Af, = Af{r)T- 

where /(r) and Af’^{r) are real, </>x,o(?') and S{r) being complex. At infinity, we have 


(4.5) 


^str 


lim f(r) = M, 
'’—>•00 

lim Af(r) = -, 

'’—>^00 g 


(4.6) 
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i.e. the value of / ensures Eq. (2.11) holds, while the gauge field cancels D^'S. After generalizing 
the notation of Eq. (2.12), the boundary conditions for the other fields satisfy 


^ Ifj^ (^) — 4^x,a,(+)^ 
""ihr Sir) = S(+). 

r^oo ' ' 


(4.7) 


At the center of the string, 


/(O) = 0, and Af"(0) = 0, 


(4.8) 


while cylindrical symmetry imposes 


—^^(0) = 0, and “ 3 ~( 0 ) = 0- 
dr dr 


(4.9) 


4.2 Lagrangian and equations of motion 

With the ansatz of the previous section, we can simplify the model tremendously. The kinetic 
term yields (with implicit summations on the representations singlet of the SM) 


K = -2 


(V^ - igAf)a - 


- - |(V^5)|2 - (4.10) 


where E®*’’ = We label (7, d) the restricted representation complex conjugate 

to {i,a). No scalar products between vectors in representation spaces are present due to the 
normalization choice of Sec. 2.3. There is no cross-terms since singlet quadratic terms can only 
be built from products of two conjugate representations. 

Writing down the potential is a bit trickier, as it contains high order terms whose derivation 
w.r.t. the fields need to be done with care. We add a subscript to the VEV indicating in which 
representation is the product we consider, {XY)^ denoting the field in the representation of Z 
coming from the contraction between X and Y. 

In the SO(IO) case, the 126 and 126 representations are fifth-rank anti-symmetric respec¬ 
tively self-dual^ and anti-self-dual tensors Yijkim and Yijkim, while the 210 is a fourth-rank anti¬ 
symmetric tensor ^ijki- The singlet which can be formed with these fields is Yijkim^ijkno^imno- 
Differentiating with respect to S, we obtain in the same representation as S, which 

is 5 [^[ij\ap^ap\kim] + i<Aijkimabcde^abap^ai3cde) [58, 59], totally antisymmetric and self-dual as 
expected. 

With these notations and with implicit summation on all the indices but e.g. I when 
considering the scalar potential becomes (noting that there is no /3-coupling, as explained in 
Sec. 3.4) 

Es = (4.11) 

for the inflaton part, 


Es = -h K^ss*aa* -h + mj:Kaa*S* + h.c. 

gCr(T*(/)a;,« + h.c. r]a;K{T.^x,a)f:^Q{^)lcra*S*4>x,a + h.c., (4.12) 

and 

Es = Es(S^S), (4.13) 

‘‘Self-duality being here defined by = jif^ijkimabcde'^abcde- 
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for the string forming field part, and 


V^j — \mij\ I (^x,q:^j/,/3)4,^^q| 4‘x,a4’x,a4^y,^4^y,f5 

+ + h.c. + + h.c. 

+ + h.c., (4.14) 

Vi!>. =V<s>j{I i — >i,Ji — >j), (4.15) 

Vij = V^j{I i —;■ i, J i —>■ j), (4.16) 

for the other fields. An example of such potential with the conventions used in this paper can be 
found in Ref. [58]. 

The Lagrangian is 

C = K-Vs-V^-V^- Y. (4-17) 

with the kinetic term given in Eq. (4.10). Note that we recover the abelian Higgs model with no 
additional numerical factor in the limit described in Sec. 3.3, due to the conventions used. 
Finally, using the ansatz given in Eq. (4.5), we obtain the following equations of motion 


'P'ia + 




±11 I 

'Pi.a + — 


dV 

dV 

d^t^la ’ 
dV 

r ’ 

r, dS*’ 

Q"-- = 2g^fQ, 

r 


2df 


^'la + 


</>U 


(4.18) 


where a prime means a derivative with respect to the radial coordinate {' = d/dr). We also 
introduced the field 

Q{r) = n- gAf^{r), (4.19) 

which is a real field function of r only with for boundary conditions 


Q(0) = n, and lim Q{r) = 0. 


(4.20) 


Eq. (4.18) with associated boundary conditions can only be solved once the actual theory 
is implemented, thus giving the relevant and necessary coefficients in the potential. 


4.3 Modification of the energy per unit length 

We now evaluate the influence of the extra fields contribution to the energy per unit length, when 
a perturbative study of the condensation of these extra fields in the core of string is possible. 

We take into account the possibility to work with high dimensional representations, of char¬ 
acteristic dimension N, as is often the case in GUTs. Let us remind that we chose the vectors 
which define the VEV directions to be normalized, i.e. (^x,o)o(^a;,o)o = 1 and (S)o(S)o = 1. 
The cubic and quartic contraction between these VEV directions must be smaller than 1 due to 
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the Cauchy-Schwarz inequality. They can be approximately estimated^ to be of order 1/y/N for 
the cubic terms, for example {^x,a^y,i3)^ order 1 /A^ for the quartic terms, as 

Important remarks must be done at this step. On the one hand, we consider here a dynamical 
case, and not only a static and uniform conhguration, where it is sufficient to independently 
ensure that the different T-terms vanish (e.g. when studying the SSB scheme). Then, the 
different multiplicities of the T-term components must be taken into account. On the other 
hand, the choice of conventions in the definition of the superpotential and the kinetic part of 
the Lagrangian (where we included here no multiplicative factor) and of normalizations for the 
constant vectors in the representation space (see Sec. 2.4) can considerably affect the formulation 
of the model. 

A complete model using the conventions of the present paper is described in Ref. [58]. In this 
SO(IO) case, taking the VEV directions for the singlets of the SM as given in Ref. [3, 5, 8 ], and 
after normalization, we hnd cubic coefficients like l/( 10 \/ 2 ) and l/( 6 \/ 6 ), and quartic coefficients 
like 1/54 and 1/164. This is roughly in agreement with what was expected with a characteristic 
dimension of order 100 for the representation, and sufficient to do a first approximation. Note 
that Ref. [58] also shows that the different formulations are indeed equivalent. 

Let us evaluate a rough order of magnitude for the potential, and for that purpose consider 
a generic held 0, without specifying indices. From Eq. (4.14), we have 


V 


2,2 , 

m cp -\ -——h 


N 


N 




VN y/N 


N 


(4.21) 


To evaluate the contributions of (p to the energy, we consider its characteristic scale of variation 
due to the presence of the string. For this purpose, we estimate the value of p at the center of the 
string and at inhnity by taking the values which minimize the potential for a = 0 and a = M, 
that we write respectively po and po + 4>i- We obtain po ~ (y/Nm)/X, and pi {r]M‘^)/{y/Nm) 
for the perturbation parameter, leading to 


pi XrjM"^ 

- , 

po Nm? 


(4.22) 


Since N ^ 1, M < m, and one has A and 77 < 1 in most of the cases, a hrst approach by 
considering the modihcations of </> as a perturbation is often relevant. 

When pi </> 0 ) we can evaluate the maximal modihcation of the energy per unit length due 
to the condensation of p in the string. Close to the conhguration with pQ, which is a minimum 
of the potential, the additional potential term is 


5U ~ 


rdr mPp'l 


Nk‘^ ’ 


(4.23) 


which hnally gives, since Uq rsj m 2 , 

6U rf 
Co "" ' n ^' 


(4.24) 


This result gives a criterion to estimate if the toy model description of the cosmic strings is relevant 
from a macroscopic point of view. Note that when pi can not be treated as a perturbation or 
when 6U <C C is not verihed, the estimate (4.23) of 5U is meaningless, and a complete calculation 
must be done. 


®For this purpose, we describe the VEV directions by vectors of N components of value 1/y/N (in order to be 
normalized), and the vectors formed from two different VEV directions by N component of values 1/N. A scalar 
product of two vectors giving N times the product of their components, we found the results used. 
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Up to now, we did not consider the inflaton S. A careful examination however shows that 
we can estimate the modifications of S due to presence of the string to be of order 4>i, and that 
the contribution to the energy per unit length of this held is of the same order or lower than the 
contribution of 4>i. It can be understood by the fact that S has no characteristic scale, which is 
necessary for it in order to play the role of the inhaton. 

5 Many-field strings 

5.1 Ansatz and boundary conditions 

We now turn to the case of a many-held string assuming the minimal structure ansatz, i.e. with 
only the restricted representations singlet under the SM taking non vanishing values. We denote 
by ^x,a the sub-representations charged under U(l)str) using as before the notations of Sec. 2.4 
for the different restricted representations. In the case of a many-held string, such representations 
appear in the GUT held content, and the superpotential contains some /3-couplings. As U(l)str 
is not broken before the last step of SSB, see Sec. 2.3, these particular helds must have vanishing 
VEVs before the end of inhation, yielding {^x,a) = 0. 

An example of such a held can be the 16 representation of SO(10), with a /3 coupling 
through 16x16x126. Its SM singlet is contained in its (1,2,4) restricted representation under 
the Pati-Salam group [SU(2)ixSU(2)/{xSU(4)]. Then, the /3-coupling can be constructed from 
this previous sub-representation and the (1,3,10) representation under Pati-Salam contained 
in the 126 representation of SO(IO). Decomposing SU(4) to SU(3)c'xU(1)s_l, the restricted 
representation of this additional held is charged under the IJ{1) b-l which can play the role of 
U(l)str [19, 59]. A similar coupling can be formed with the 16 and 126 representations. 

More than one held and its conjugate are now charged under U(l)str, so we have to specify 
these charges. We call them qx^a and gs = —g^. As we can multiply these charges (which can be 
fractional) by a common factor without loss of generality, we choose the smallest set of charges 
where they all are integers, i.e where a rotation of 2tt in the U(l) group is the smallest one which 
reduces to identity. It will be convenient in what follows to dehne the winding number. 

The D-term condition given in Eq. (4.3) should now include all the helds that are charged 
under U(l)str- This condition is a priori not sufficient to ensure that a* = a and similar relations 
for all the helds charged under this abelian symmetry, i.e. (i>*x a = 4>x,a, with ^x,a in the conjugate 
representation of ^x,a- We will however assume it from now on, since this is a solution of Dstr = 0, 
and so a global minimum of the associated potential. 

The same topological arguments than above imply that all the helds which are charged 
under U(l)str vanish at the center of the string [II] : 

{ (S)(^=o) = (^)(r=0) =0) /Cl') 

\ (^x,a)(.=0) = 0. 

We consider the following ansatz for an abelian cosmic string associated with the generator 

j-str 

a = Ury^^^\ 

^x,a = 

— 4^x,a{'^')j (b'2) 

S = S{r), 

A, = Af{r)T^^^5f,. 

In these equations, the fx and Af"’^ are real functions, and 9x,a are reals constants. The possibility 
to freely dehne an origin for the coordinate 6 permits to put no initial phase to the held a. 
The integer n is well identihed with the winding number, taking into account the choice of 
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normalization we did for the charges. We also have the following boundary conditions, similarly 
to the single-held string case, 


lim f(r) = M, 

>-CXD 

Jil^ (f>x,a (^) — 4’x,a,(+)j 

lim ,S(r) = ,S(+), 
lim Af(r) = -, 

"^oo ^ ^ g ' 


(5.3) 


at inhnity, and 


/(O) = 0, = 0, 

^(0) = 0, and 

dr 


( 0 ) = 0 , 


(5.4) 


at the center of the string. 

At inhnity, our ansatz simply gives an absolute minimum for the potential on which we 
applied a local U(l)str transformation. It is the standard restated result for a cosmic strings, but 
here with several helds charged under the string forming U(l). It shows that the angles Ox^a are 
not freely chosen, but are those which give the absolute minimum of the potential and thus dehne 
the SM symmetry at inhnity, in addition to the limits given in Eq. (5.3). 


5.2 Modification of the energy per unit length 

We now turn to the modification of the energy per unit length from standard toy models due to 
the condensation of the additional fields in the core of the string, with a perturbative approach. 
As in the case of a single-field string, we work with a generic field 4>, without considering its 
indices anymore. 

The potential term gives 


y + 


w 


mr](l)a‘^ mfiacfP' 

~ir ^ ^ ^ 


y/N y/N 


^ m\<p^ ^ ^ ^ /3Xa(p^ 


y/N 


N 


N 


N 


(5.5) 


Introducing ~ {y/Nm)/X, two perturbation scales appears, (j)i ~ {r]M^)/{y/Nm) and (j)'i ~ 
{(5M)/X. So, in order to make a perturbative study of the string, the small parameters we have 
to consider are 


(h 

4>o 


XrjM'^ 


and — 


4>o 


PM 


(5.6) 


Nm"^ 4>0 y/Nm' 

If a perturbative study is possible, we then estimate the maximal modification of U due to the 
second term to be of order 


6U' ~ J rdr rr/(f)\ 




(5.7) 


which gives 

6U' 

We see that even in a model where the coupling constants are smaller than 1, we cannot in general 
consider that the modification of the energy per unit length of the strings can be treated as a 
perturbation. It means that in the case of a many-field string, it is necessary to do a complete 
study of the microscopic structure of the string, and also that the toy model approximation is in 
general not valid. 

Finally, we see that the simplification of Sec. 2.1 consisting in omitting terms like 
and can be justified a posteriori. Indeed, these terms must be studied in the case of a 
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many-field string, as the field $ appearing in this cubic term has to be charged under U(l)str 
in order to provide a singlet term. However, its contribution to the energy will be similar to a 
standard term, which is already taken into account here. 

6 Conclusions and discussions 

Let us summarize by stating the different results obtained in this paper, focusing on the abelian 
strings forming at the last step of the SSB scheme, at the end of inflation. We worked in the 
framework of a SUSY GUT with hybrid inflation, but some results are more general, and we state 
for each point which hypotheses are actually necessary. 

i) The minimal structure model for realistic cosmic strings contains at least all the fields taking 
non vanishing VEVs at the end of the SSB down to the SM, which all condense in the core 
of the strings. These fields are singlets of the SM. This property neither depends on the 
GUT nor on the inflationary process. 

ii) It is possible to distinguish two classes of strings, the single-field strings where only one 
field and its conjugate are charged under U(l)str) and the many-field strings. The way the 
field content performs the SSB scheme is sufficient to determine which kind of strings form, 
independently of having a SUSY model or of the inflationary process. Different couplings 
appear in each model, which may give different phenomenologies for the strings. 

iii) For each class of strings and with a SUSY GUT and hybrid inflation, we gave an ansatz 
and the associated boundary conditions to describe the minimal structure of the strings, 
completely defining it from a mathematical point of view. These strings are singlet of the 
SM. The fields take intermediate values in the core of the string between the minimum of 
the potential defining the SM at infinity and the configuration taken by the fields before 
the end of inflation. This ansatz can easily be generalized for other models. 

iv) Going from the GUT description to the minimal structure ansatz, important numerical 
factors appear and must be taken into account. How these factors appear depends on the 
conventions and normalizations used. A special care has to be taken in SUSY models, 
where these factors are omitted when studying the static configurations which minimize the 
potential, i.e. with all F-terms independently vanishing. We emphasized how to normalize 
the GUT in order to recover the abelian Higgs model in the limit where the couplings 
between the string-forming Higgs and the other fields go to zero. 

v) In the case of a SUSY GUT with hybrid inflation, we performed perturbative estimates of 
the modification of the energy per unit length with respect to standard toy models, which 
are given in Eqs. (4.24) and (5.8). These results are very different in the case of a single¬ 
field string, where the modifications are sizable in a high-coupling limit, and in the case of 
a many-field string, where the modifications are always important and require a complete 
computation of each model. 

These results represent a first step of a more thorough investigation of cosmic strings taking 
into account their realistic structure. A complete study of the microscopic structure and the 
energy per unit length of such strings in a given SO(IO) model has been performed in parallel, 
including numerical solutions, and can be found in Ref. [58]. As discussed in Sec. 3.2, several 
properties of the strings can also be modified in this framework and have major cosmological 
consequences. Indeed, the Higgs fields condensing in the core of the string could carry bosonic 
currents [33-37]. Moreover, their superpartner could carry fermionic currents through their zero 
modes [33, 38-42]. Furthermore, this complex microscopic structure could qualitatively modify 
the intercommutation process [43-48] and thus the evolution of the cosmological string network. 
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hence modifying its consequences on the CMB [10, 61-65]. The cusps evaporation [49-53] also 
deserves more thoughts, since in our case, knowledge of the fields present in the string core implies 
knowledge of the relevant branching ratios into specific particles. 
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